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I.  Executive  Summary 

A  multiple-element  scale/coherence  decomposition  of  the  Navier-Stokes  equations  is 
employed  to  compute  the  Lagrangian  evolutionary  dynamics  of  the  dominant  three- 
dimensional,  non-linear  vortical  features  in  high  Reynolds  number  wall-bounded  turbulent 
shear  flows.  Organized  turbulent  and  wave-like  structures  of  various  scales,  their  inter¬ 
actions  and  transitions  between  them  are  computed  distinctly  and  subsequently  shall  be 
modeled.  The  resultant  moael  is  denoted  'coherent  structure-reflective'  since  it  is 
derived  from  equations  whose  prediction  of  organized  turbulence  elements  (evolution/ 
propagation/occurrence)  will  have  been  validated  through  comparison  to  observations. 

This  analysis  is  carried  out  initially  for  constant  pressure,  two-dimensional  boundary 
layers.  The  complex  effects  of  a  solid  boundary  (and  its  associated  low  Reynolds 
number  flow  regime)  on  the  difficult  instability  of  a  non-inflectional  mean  turbulent 
velocity  profile  indicate  the  virtues  of  distinguishing  between  elements  of  various 
scale  and  coherence.  In  addition,  the  work  points  out  the  manner  in  which  this  approach 
facilitates  treatment  of  all  the  original  N-S  equation  non-linearities  via: 

1.  The  decomposition  itself, 

2.  Perturbations  analysis,  and 

3.  Small  scale  stress/flux  modeling. 


Subsequent  work  on  complex  shear  flows  involving: 


% 


1.  Streamwise  and  normal  pressure 
gradients  due  to  curvature, 

2.  Unsteady  external  flow, 

3.  Three-dimensionality,  and 

4.  Compressibility, 

is  also  amenable  to  structural  computation  and  modeling  which  is  in  accord  with 
observation  and  the  following  schematic. 
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II.  Introduction 


The  accurate  prediction  of  untested*,  high  Reynolds  number,  wall-bounded,  turbulent 
shear  flows  in  complex  situations  is  sought  from  a  new  class  of  representations  of  the 
time-  or  ensemble-averaged  flow.  Such  theoretical  formulations  are  herein  denoted  as: 


"Coherent  Structure-Refiective" 


in  order  to  emphasize  the  fact  that  the  various  Eulerian  closure  formulations  embody 
current  experimental  reality  and  are  derivable  from  equations/computations  of  the 
time-dependent  organized  turbulent  structure  and  its  Lagrangian  evolutionary  details. 

Note  that  such  analyses  are  important  in  their  own  right,  apart  from  subsequent  model¬ 
ing,  in  the  understanding  of  non-linear  vortical  dynamics.  In  fact,  the  simplest  set  of 
equations  which  contain  the  essential  structural  properties  required  to  predict  turbulent 
transoort  are  tested  by  direct  comparison  with  observations.  Following  verification, 
the  governing  structure  equations  shall  be  processed  to  provide  a  model  which  resides 
somewhere  between  current  multi-equation  modeling  formulations  and  LES/full  N-S 
equation  solvers  (which  provide  'typical'  turbulent  realizations  that  may  be  ensembled), 
benefiting  from  developments  in  both  fields. 

For  the  foreseeable  future,  at  flight  Reynolds  numbers  of  interest,  modeling  will  be 
required.  However,  it  is  the  intent  of  this  work  to  develop  formulations  more  nearly 
in  accord  with  experimental  observations.  Improvement  in  the  prediction  of  these 
fiowfields  is  acnieved  by  addressing: 

♦To  the  extent  that  insufficient  data  is  available  for  traditional  case-by-case  turbulence 
model  'tuning'. 
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Multiple-scale, 

Anisotropy,  and 
Non-gradient  transport, 

phenomena  directly  in  the  models,  in  a  manner  compatible  with  experimental  data  and 
rigorously  derived  from  the  deterministic  structure  dynamics  equations  developed  here. 
In  addition  (and  particularly  for  the  complex  flows  to  be  addressed  in  subsequent  work), 
the  approach  is  extended  to  accomodate  truly  coherent  wave-like  oscillations  due  to 
curvature,  unsteady  excernal  flow,  three-dimensionality,  compressibility  and  low  Rey¬ 
nolds  number  (near-wall)  effects. 

As  a  result,  the  overall  procedure  is  a: 

"Multiple-element  Turbulence  Decomposition", 

where  uniquenesses  of  the  flowfield  structures  associated  with  various  scales  and 
coherences,  their  interactions  and  transitions  between  them  are  computed  distinctly 
(and  subseauently  modeled)  according  to  the  following  prescription: 


STEADY  N SCALES 

PERIODIC 


INTERACTIONS _  M  SCALES 

TRANSITIONS  INCOHERENT 
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In  this  methodology,  all  nonlinearities  are  accomodated  by: 


1.  The  decomposition  itself, 

2.  Perturbation  methods,  and 

3.  Modeling  small  scale  stresses/fluxes, 

resulting  in  an  expanded  set  of  weakly  non-linear  equations  whose  soultion  structure 
is  more  readily  accessible  than  that  of  the  primitive  or  sub-grid-scale-modeled  large- 
eddy  equations. 

The  role  which  observational  data  or  LES/full  N-S  solutions  play  in  calibrating  or 
verifying  this  approach  is  more  subtle  than  in  traditional  turbulence  modeling.  The 
primary  difficulty  in  this  regard,  even  for  two-dimensional  constant  pressure  boundary 
layers,  is  that  conventional  modeling  is  rooted  in  fundamentals  which  are  not  in  accord 
with  observation.  As  a  result,  extensive  empirical  adjustment  of  the  various  'constants' 
and  functionals  is  required  to  post-diet  the  experiment  from  which  the  models  were 
'tuned'.  The  need  for  calibration  is  not,  however,  the  Achilles'  heel  of  modeling.  The 
problem  is  the  situ  tion-specific  tuning  of  models  dictated  by  the  faulty  premises  on 
whicn  they  are  based.  Sufficiently  general  models  may  require  'tuning'  but  oniy  in  a 
sing.-.-,  simplified  flow;  ala'  the  hot-wire  anemometer  calibration. 

Por  this,  reason,  it  is  eaually  as  difficult  to  stretch  any  single-fluctoating-element 
two-equation  (e.g.  k-e)  model  into  being  'coherent  structure-reflective'  by  comparing 
it  with  foe  mijitiole-element  computational  results  presented  here.  This  is  not  surpri¬ 
sing  since  t-ese  models  are  grounded  in: 


1.  Essentially  isotropic, 

2.  Single  scale,  and 

3.  Gradient  transport 

notions.  (In  contrast,  free-mixing  situations,  excepting  low  Reynolds  number  regions, 
are  essentially  single-scale  for  which  t-ou'H.  9u/3y  is  not  a  gradient  transport 
assumption.  In  these  flows,  differential  equation  representation  of  the  length  and 
velocity  scales  enhances  already  good  predictions  by  accounting  for  non-equilibrium 
effects.!  In  addition,  current  Reynolds  Stress  Tensor  models  are  not  sufficently  mature 
for  reliable  application  to  complex  flows  or  adaptation  to  coherent  structure  reflecive 
status.  Tnerefore,  the  bulk  of  the  work  reported  here  focuses  not  on  adjusting  other 
models  but  on  computing  coherent  structures  and  formulating  models  appropriate 
thereto.  The  associated  validation  of  coherent  structure  reflective  models  is  a  simple 
non-iterative  comparison  with  observation  or  LES/N-S  solver  results,  once  it  has  bee 
established  that  the  deterministic  structural  equations  are  grasping  the  essence  of  the 
three-dimensional,  non-linear  vortical  dynamics.  In  this  regard,  some  caution  must  be 
exercised  when  comparing  computed  streakiines  to  observations  and  then  making 
inferences  about  the  _agrangian  structures  involved. 

Apart  from  the  discrepancy  between  single-scale  modeling  concepts  and  multiple-scale 
fiows,  another  phenomenon  arises  in  wall-bounded  situations  which  adds  justification 
to  the  scale/coherence  decomposition  approach  developed  here.  This  is  associated  with 
the  low  Remolds  number  nature  of  the  near-wall  region.  In  fact,  the  structure  of  this 
unsteady  viscous  flow  Is  more  nearly  instability  wave-like.  Therefore,  when  representing 
ail  appropriate  f!uetuations  in  the  flow  near  a  wail,  both  wave-like  and  turbulent 


□recesses  may  be  required.  The  solid  boundary  essentially  damps  the  turbulent  fluctua¬ 
tions  wnile  making  tne  organized  oscillations  highly  anisotropic,  clearly  a  difficult 
nature  for  essentially  isotropic  modeling  to  accomodate  under  the  single-element  mo¬ 
deling  umbrella.  In  fact,  this  is  the  reason  that  e  is  set  equal  to  zero  at  the  wall, 
wail,  whereas  the  total  dissipation  is  finite  due  to  anisotropic  effects.  Therefore,  even 
ignoring  the  multiple-scale  nature,  conventional  single-element  modeling  oversimplifies 
the  near  wail  Huctuaticns/interactions/transitions.  Similar  complications  arise  at  the 
fluctuating  interface  with  the  outer  flow. 

These  phenomena  indicate  tne  need  to  apply  multi-element  modeling  to  even  simple 
flows*.  However,  the  two-dimensional,  constant  pressure  turbulent  boundary  layer  is 
an  appropriate  starting  point  for  the  deterministic  structural  computations  of  interest 
here.  This  is  because  tne  turbulent  mean  velocity  profile  is  non-inf lectional  and  thought 
to  oe  stable,  even  to  'aooarent  viscous'  instabilities  found  in  laminar  flows.  It  will  be 
shewn  that  the  resolution  of  this  problem  is  associated  with  a  delicate  interplay  between. 

1.  Mean  velocity  profile  details, 

2.  Small  scale  stress, 

a.  non-equilibrium 

b.  anisotropy 

?.  ’Son-parallelism,  and 
4.  Mon-linearity, 

excused  o>  tKe  stale  ss^ere'ce  decomposition. 

*Tne  final  sir*  Remolds  ^u^ber)  stages  of  even  sheared  grid  turbulence  develops 
highly  anisotropic,  wave-like  structures. 


Iv.b.  The  Prescription 

•n  croer  t:  implement  the  pnvsics  of  Section  iv.a.  into  a  comprehensive  mathematical 
-Tiociei,  the  ‘cliowinc  prescription  has  been  devised. 

i.  Decomocsition. 


Multi-element  decompose  the  equations  of  motion  to  differentiate 
between  events  of  various  scale  and  coherence,  their  rter act  ions 
and  transitions  between  them.  Accomodation  for,  and  subsequent 
modeling  of  the  orimitive  non-linearities  is  therely  faciiitatec 
through; 

a.  The  decomposition  itself, 

fc.  Perturbation  expansions,  and 

c.  Modeling  the  small  scale  stresses./fluxes, 

at  the  orice  of  an  expanded  (yet  simplified  and  essentially  rigorous) 
set  cf  governing  equations  as  follows: 


Complex  flow. 


Wave-turbulence  interactions  and  transitions.  Large  scale- 
small  scale.  Signal-background  6-6^  (non-parallelism)  Remnant 
'summing'  leads  to  various  added  elements  u'Q  and  u'g  ,  etc. 

Self-perpetuation. 

Remnant  summing. 
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Interactions 

1.  The  superlayer. 

High  speed  (U)  external  flow  'runs  into'  high  V  ejection  jets  propa¬ 
gating  at  U/H,  partially-entrained  and  partially  disturbed,  possibly 
setting-up  the  downstream  flow  with  precursor  information.  Upstream 
and  downstream  'radiation'  occurs.  The  interface  between  the  distur¬ 
bed  external  irrotational  flow  and  the  rotational  shear  flow  is  distinct 
and  a  clear  transition  between  differing  fluctuating  states  (coherent 
vs.  incoherent?).  The  time-average  view  of  a  smoothly-varying 
vorticity  is  irrelevant. 

2.  The  viscous  sublayer. 

The  region  is  unsteady  and  low  Reynolds  number,  indication  a 
distinctly  non-chaotic  nature.  It  may  be  understood  in  terms  of 
two-layer  (i.e.  spatially-varying  v)  Langmuir  roils  and  a  streaming 
phenomena  as  a  result  of  wave-turbulence  interactions. 

3.  Inner  flow-outer  *"!ow  synchronization. 

"he  c:°ssure  field  and  the  external  flow,  which  passes  over 
Durst  events  anp  proceeds  downstream,  may  serve  to  unify  the 


u  «  f  •  T*  7  *  U  »  !■ 


Vorticity 


1.  Three-dimensional  vortex  tilting/stretching  and  breakdown. 


2.  Vertical  vorticity  -  vertical  velocity  resonance/non-linear  dispersive  wave 
mechanics. 


3.  Three-dimensional  local  separations. 


Curved  and  sheared  flow  around  'jets'  associated  with  the  fundamental 
instability  may  introduce  another  class  of  hairpin  vortex  structure. 


4.  Langmuir  circulation  patterns  in  flows  with  spatially  non-homogeneous  v. 


3.  Goertler  vortices. 


hree-dimensionality 


1.  Three-dimensional  boundary  layer  instability. 


2.  Three-dimensional  local  separations. 


3.  vertical  vorticity  -  vertical  velocity  resonance. 


u.  Heiicity. 


5.  Self-perpe‘.uation. 
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2.  Critical  level  interactions  are  significant  to  understand  the  fundamental  instability. 

'Viscous'  and/or  non-linear  effects  are  important  but  the  initial 
value  problem  nature  of  the  vortex  dynamics  evolution  may  also 
impact  the  situation. 

3.  Multiple  scale  perturbation. 

The  fundamental  instability  drives  a  free  shear  layer  inviscid 
instability  involving  large-scale/small-scaie  and  coherent/ 
incoherent  interactions.  There  is  an  equivalence  between  large- 
scale  (residual)  non-linearity  effects  and  the  effect  of:  coherent 
fluctuations  on  the  large  scale,  large  scale  fluctuations  on  the 
small  scale,  etc.,  leading  to  "multiple-scale  perturbation" 
amplitude  envelope  evolution  problems. 

A.  Cause/effect  vs.  non-linear  flowfield  unity. 

Note  tnat  sequential  cause-effect  within  a  particular  scale/ 
coherence  is  equivalent  to  the  aforementioned  effect  of  coherent 
fluctuations  on  large  scales,  large  scales  on  small  scales,  etc. 

5.  Self-perpetuation. 

In  a  linear  system,  no  Qualitative  change  in  dynamics  would  occur 
as  a  result  of  turbulent  burst  structure  remnants  combining. 


carry  their  own  noise  generator.)  However,  like  laminar  flow,  in 
terms  of  the  given  local  state  (e.g.  6),  the  typical  details  can  all 
be  specified.  But  determining  the  local  state  requires  initial  condi- 
tior  and  boundary  condition  (irradiation)  statistics  and  history, 
including  streamwise  inhomogeneities  and  non-equilibrium.  Note 
also  that  upstream  (initial  condition  contamination)  and  down- 
streami  (precursor)  radiation  of  information  is  possible  due  to  the 
fact  that  the  external  disturbed  (yet  not  entrained)  flow  is  traveling 
at  U  -  u/H  relative  to  the  L  structure. 

Non-linearity 

1.  Small  scale  non-linear  terms  are  modelable. 

According  to  Landahl,  the  large-scale  structures  are  weakly  non¬ 
linear  as  long  as  the  (background)  small-scale  'apparent'  stresses/ 
fluxes  or  sources/sinks  are  accounted  for  in  the  various  interactions/ 
transitions.  Very  fine  scale  effects  may  simply  reduce  the  effective 
Reynolds  number  but  caution  must  be  exercised  when  adding  vt  to  v. 
Clearly,  the  small-scale  non-linearities  are  non-equilibrium,  possibly 
anisotropic,  not  necessarily  Newtonian  and  Reynolds  number- 
dependent.  vf  gradients  (ut  and/or  c)  lead  to  streamwise  vorticity, 
reminiscent  of  _arigmuir  circulations. 


y  - 


The  initial  condition  variability  and  selective  amplification  of 
ambient  or  seif-induced  radiation  also  have  counterparts  in  the 
rain-field.  The  abstract  result  is  a  non-linear  box  with  multiple 
feedback,  pure  tone  excitation  and  expansion/compression  (i.e. 
amplitude-dependent  amplification). 

3.  Propagation  dynamics  of  large-scale  structures. 

The  streamwise  propagation  of  large-scale  structures  is 
understandable  in  terms  of  the  equation: 

L,t  +  (U/H)L,x  =  fL  +  g  , 

which  is  the  familiar  time-dependent  integral  TBL  equation  if  L  =  0 
and  g  =  Cf  .  The  resultant  space-time  focusing  and  wave-like 
i.vs.  diffusive)  nature  of  these  non-linear  dispersive  phenomena  is 
accessible  from  such  an  equation. 

4.  Observed  'complex'  st:  ;3klines/pathjines/vortex  lines  may  arise  from  simple 
streamlines,  thus  pointing  out  the  importance  of  the  right  viewing  reference. 

3.  Independence  of  initial  concitions. 

If  the  mitia.  conditions  (for  the  wave  instability  and  coalescence) 
are  the  cm/  source  of  irregularity,  the  solution  is  not  strictly 
independent  oc  initial  conditions.  (Note  that  all  numerical  schemes 


'viscous'*  and,  possibly,  non-parallel  and  non-linear  (threshold/ 
hypteresis)  effects.  Note  that  experimental  data  may  suffer  some 
spatial/temporal  'smearing'  due  to  the  signal  processing.  Subsequent 
to  the  discovery  of  a  fundamental  mean  profile  distortion  instability 
is  a  secondary  free-shear  layer  inviscid  instability,  due  to  non-linear 
mean  profile  distortion,  resulting  in  a  multiple-scale  perturbation 
amplitude  evolution  problem. 

2.  Chaotic/statistical  occurence  of  organized  structures. 

Ala'  the  rain-field  model,  this  may  be  the  natural  mathematical 
state  of  weakly-constrained  continuous  vector  fieids,  as  indicated 
by  strange  attractor  theory.  Clearly,  our  Eulerian  point  of  view 
is  deceiving  and  unnecessarily  complicating,  as  the  raindrop  model 
makes  clear.  (An  alternative  metaphorical  view  of  turbulence  is 
related  to  traffic  flow  which  points  up  the  need  for  a  Lagrangian 
approach.  Conditionally-sampled  data  is  inevitably  an  average  over 
various  typical  turbulent  structures  and  distinctions  between,  e.g., 
trucks,  buses  and  cars  are  smeared  out.)  Turbulent  flow  'patchiness' 
is  clearly  more  readily  described  by  fractals  or  pdf’s  which  depend 
on  (flow)  historical  inputs  regarding  instability  thresholds  and 
hysteresis,  size/speed  variability  and  a  metastable  nature. 
Interactions  between  burst  structure  remnants  into  a  disturbance 
o*  Sufficient  ’size’  is  represented  by  the  various  non-linear  terms. 

•Note  gradients  in  u  and/or  d  are  secondary  vorticity-inducing  and 
stream* ise  vorticity  acts  like  enhanced  transport  to  the  large 
scale  'low. 


I 

I 

This  process  continues  as  long  as  the  raindrops  have  potential  energy,  and  there  is  j 

no  surprise  that  perpetuation  and  statistically-scattered  impact  occur.  Similarly,  ! 

once  turbulent  flow  instability  waves  coalesce  (irregularly,  dependent  on  statistically-  ! 

t 

distributed  nucleation  sites),  their  propagation/interaction  and  bursting  instability  ^ 

proceeds  as  long  as  the  main  flow  has  kinetic  energy  and  is  unstable. 

Followinc  the  rain-field  analogy  a  step  further,  it  is  sufficient  (for  the  composition  of  J 

a  comprehensive  model)  to  assess  deterministic  raindrop  propagation,  interactions  and 
instability  following  a  given  organized  structure.  Subsequently,  it  is  possible  to  assemble 
a  statistical  model  of  drop  and  burst  occurence  based  on  variability  in  the  coalescence  ^ 

initial  conditions  and  various  threshold/hysteresis  effects  in  the  individual  droplet  dynamics. 

Returning  now  to  the  non-global  details,  several  features  nf  the  TBl  flow  are  enumera¬ 
ted  here,  in  accord  with  the  previous  definition,  which  are  part  of  the  long-range  focus 
of  this  research: 

Time-dependence 

1.  Periodically-infiectional  conditionally-sampled  mean  TBL  velocity  profile. 

The  time-averaged  mean  profile  is  non-infiectional  and  thought  to 
be  stable.  Therefore,  primary  instability  must  be  related  to  the 
small-scale  non-linear  terms  (stresses/fluxes  or  sources/sinks)  ala' 
laminar  boundary  layers.  Conditionally-sampled  structure  data  should 
be  compatible  with  computations  utilizing  non-equilibrium/non-uniform 
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RAIN-FIELD  'TURBULENCE' 


3- 


The  Building  Blocks 


When  attempting  to  assemble  a  comprehensive  representation  of  turbulence  which 
is  in  accorc  with,  and  reflects  organized  structure  physics,  it  is  appropriate  to 
catalog  and  examine  tne  various  proven  or  proposed  component  processes.  It  may 
seem  trite,  but  it  is  useful  for  organization,  to  define  turbulence  as  complex  time 
Dependent,  non-linear,  three-dimensional,  vortical  interactions.  In  the  process  of 
documenting  the  content  of  the  multi-element  scale/coherence  approach,  various 
phenomena  which  are  confronted  fit  within  one  or  more  of  these  categories. 


Overlaying  many  of  these  phenomena  is  a  global  view  of  turbulence,  similar  in  some 
regards  to  the  following  raindrop/rain-field  scenario: 


IV.  Results 


The  research  wor*  presented  here  involves  analytical/computational  investigations 
cf  various  aspects  of  the  evolution  (in  the  convected  frame),  propagation  and 
statistical  occurrence  of  large  scale  structures  in  a  TBL  flow.  Prior  to  the  presen¬ 
tations  of  calculations,  sections  which  describe  the  building  blocks  of  the  analysis, 
tne  computational  prescription  and  an  enumeration  of  the  various  facets  of  the 
mathematics  and  physics  which  this  approach  seeks  to  embrace  are  presented. 


The  secondary  objective  of  this  work  is  to  lay  the  foundation  for  multi-element  com 
putations  of  turbulent  structure  and  subsequent  modeling  of  more  complex  flows. 
Specifically,  the  intent  is  to  address  the  effects  of: 


1.  Longitudinal  curvature-  including:  strong  P  and  P  ;  Possible 

re-transition/wave^  and  localized 
separation  at  the  leading  edge;  non¬ 
boundary  layer  wake  interaction/ 
separation/Goertler  vortices  formation 
at  the  trailing  edge, 

2.  Unsteady  external  flow- impulsive,  periodic  and  stochastic, 

3.  Three-dimensionality-  including  separation  bubbles,  and 

h.  Compressibility-  vorticity/dilatation/pressure  fluctuation 

interactions  of  various  scales, 


These  phenomena  shall  be  superposed  (not  necessarily  linearly  of  simply)  on  the  two- 
dimensional  constant  pressure  boundary  layer  has  been  studied  preliminarily  in  this 
first  year  of  research. 


III.  Objectives 

The  primary  objective  of  this  multi-year  research  effort  is  to  compute  turbulent  struc¬ 
tural  characteristics  and  integrate  them  into  models  of  the  flow  such  that  accurate 
predictions  of  untested  situations  are  possible  without  case-by-case  'tuning'.  Tradition¬ 
ally,  Eulerian  flowfield  predictions  focus  directly  on  computing  surface  stress  (including 
pressure  for  viscous  interaction  cases)  and  flux  distributions,  from  which  they  were 
tuned,  such  that  the  forces,  heat  loads  and  wake  signature  characteristics  may  be  com¬ 
puted.  Ideally,  the  modeling  may  eventually  become  sufficiently  general  to  accomodate 
local  or  massive  separations  and  all  of  the  complex  effects  of  interest  in  the  subse¬ 
quent  phases  of  this  research. 

In  order  to  achieve  a  reliable  model  of  the  turbulence  which  provides  such  predictions, 
we  seek  to  compute  Lagrangian  turbulence  structural  details  (evolution/propagation/ 
occurrence),  building  the  foundation  of  the  turbulent  stress/flux  representations  with 
the  requisite: 


1.  Multiple  scale, 

2.  Anisotropic,  and 

3.  Non-gradient  transport 

nature  and  accomodating  both  turbulent  and  non-turbulent  fluctuations.  Subsequently, 
ability  to  fabricate  models  directly  from  the  multi-element  structural  equations  shall 


be  assessed. 


TYPICAL  POUR  -  ELEMENT  SCALE/COHERENCE  DECOMPOSITION 
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Solve  large  scale  structure  evolution  equations  and  compare  results  to  experi¬ 
mental  streaklines. 


3.  Compute  'details': 


a.  Primary  distortion  and  secondary  free  shear 
layer  instability  amplitude  evolution. 

b.  Sublayer  vortices  -  wavy  and  acting  as 
effective  variable  v  in  small  scale  stresses? 
Superlayer  waviness  -  entrainment. 

c.  Vertical  velocity  -  vertical  vorticity 
interaction/resonance. 

d.  Non-linear  wave-turbulence  interactions/ 
transitions. 

e.  Three-dimensional  separation-type  vorticity 
generated  by  curved  sheared  flow? 


A.  Evaluate  propaqaticna!  characteristics  of  large-scale  structures. 


3.  Evaluate  statistical  occurrence:  pdf  and  fractals. 


6.  *Average/mode!  the  verified  large-scale  structure  governing  equations  based  on 
structural  event  computations. 


7.  Compute  the  mean  flow  with  a  coherent  structure  reflective  stress  model. 


3.  Evaluate  mcde.’i'q  sensitivities. 


♦item*;  c  -  I ?  are  the  sub'ect  of  on-ooino  research. 


Iv.c.  Computed  Details 


Durinc  the  first  year  c‘  work  on  this  research  the  bulk  of  the  effort  has  been  directed 
at  f'e  various  analytical  formulations  and  numerical  codes  required  to  implement  the 
ore\ iously-discussed  cresnriptior..  The  preliminary  results  to  be  discussed  shall  be  con¬ 
siderably  supplemented  in  future  research  reports  as  the  aforementioned  research  tools 
are  re'inec  anc  made  mere  productive.  The  scope  of  the  work  beina  carried  out 


i.  _arge  Scale  Structure  Analysis/Computations 

li.  nathline/Streakline  Computations 

III.Generic  Transport  Eauation  Analysis/Computaticns 

Aif  recard  to  the  item  the  first  phase  of  the  large  scale  structure  computations 
is  complete,  utilizing  the  multi-element  sea. coherence  decomposition  approach. 

The  computation  of  large  scale  structures  has  been  carried  out  utilizing  the  following 
vorticitv  equations: 
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PA-ttalLv  -  LINEARIZED  3-D  VORT'CITV  EQUATIONS 

|{u-c)  ♦  I  [(Re_1+  e) (D2-k2)  ♦  (D2c)]|<L1«  -{(Du)  +1  (De) (2D2-k2)}u3 

f 

-  i  {(De)D  -  (2D2t)}u2 

(u-c)  ♦  [(Re  *+  e)(D2-k2)  ♦  (Dc)D]j<*>2  ■  -  £  (Dujuj 

(7-c)  ♦  I  [(Re_1+  c) (D2-k2)  ♦  (D2c)]|w3  -  {(Du)  ♦  1  (Dc)  (2D2-k2)Ju1 

♦  {(De)D  -  (2D2c)}u2 

♦  i  [D2u  ♦  (Du)D]Ju2 


*~ie r.e  eq.jaticm  are  oerived  from,  the  previously-presented  scale/coherence  decom¬ 
posed  equations  and  accomodate: 

a.  spatially-varying,  non-equilibrium  small-scale  processes 

b.  multiple-scale  perturbation  analysis  of  residual  non¬ 
linearity  and  wave-like  flow  component  distortion 
ejects  (see  discussion  on  131.) 

c.  evaluation  of  vertical  velocity-vertical  vorticity 
inter  act  ion/resonance. 

Tne  following  results  for  a  typical  mode  (including  only  small  scale  process  non-line- 
arity)  illustrate  the  amplitude  and  phase  distributions  of  the  vorticities  and  vertical 
velocity  in  the  absence  of  waves  and  large  scale  structure  non-linearity.  A  simplified 
Newtonian  model  is  used  for  n..  One  developing  hypothesis  of  the  work  is  that  phase 
coherence  (at  a  stationary  phase  value)  at  various  locations  across  the  flow  is  indi¬ 
cative  or  significant  interactions  which  the  pathlines  should  reflect.  Also,  the  modal 
assemblage  (in  a  manner  which  equates  temporal  and  local  spatial  average,  or  overall 
spatial  average  for  parallel  flow)  could  be  guided  by  such  coherence.  The  frequency, 
wave  number  and  spatial  location  in  the  shear  layer  of  the  most  unstable  mode  is  in 
agreement  wit" 


experimental  data. 


"he  second  ares  of  stuay  is  the  generation  and  interpretation  of  three-dimensional 
pst-dines  and  streaklines  from  the  results  computed  in  1.  A  typical  result  which 
presents  a  two-dimensional  cut  through  the  flow  is  presented  here.  Ideally,  com¬ 
parison  between  experimental  and  theoretical  streaklines  will  serve  to  validate  the 
analysis  and  interpret  the  data. 
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rinaiiy,  the  following  generic  transport  equation: 


L,t  -  (U/H)l,x  =  fL  +  g  , 

where  g  exhibits  a  weak  inverse  dependence  on  L,  is  being  analyzed/computed  with 
respect  to  the  role  it  plays  in: 

a.  amplitude  evolution  of  secondary  instabilities- 
residual  non-linearity  plus  wave-like  distortion, 

b.  occurrence  statistics, 

c.  propagation  dynamics,  and 

d.  general  transient  boundary  layer  behavior. 


By  utilizing  a  multiple-scale  perturbation  approach,  the  secondary  instability  amplitude 
evolution,  for  both  large-scale  non-linearity  and  the  presence  of  waves,  is  found  to 
obey  the  following  equation: 

S,_  +  (?f/c?u))-1  S,x  =  aS  +  bS* 


In  addition,  current  work  on  pdf  representations: 


sP  a  bP  .  dP  .  r,  . 

b  —  + - —  - - (>u,  - - P  —  P 

dr  dx,  0^  dx.  Ox,  rf 


*  ~  v5  [ie &z.Pflz' )  )PU'.Z 


leads  tc  a  similar  formulation.  Finally,  if  L  is  identified  with  any  boundary  layer  inte¬ 
gral  thickness,  the  equation  may  be  used  to  investigate  the  propagation  of  ’bulges’  in 
the  flow  and,  subsequently,  3-D  interactions  between  laterally-displaced  structures  which 
are  advecting  principally  in  the  streamwise  direction,  /it  present,  computations  are 
focusing  or,  obtaining  stable  solutions  with  large-amplitude  disturbances  for  which 
traditional  stability  theory  is  not  rigorously  valid. 


V.  Conclusion 


On-going  work  involves  extensive  exercise  of  the  tools  described  here  and  the  inclusion 
cf  the  presence  o*  wave-like  elements  and  large-scale  non-linearity  through  a  multipie- 
scale  perturbation  approach.  In  addition,  the  optimal  representation  of  r. .  is  being 
studied  in  preparation  for  follow-on  research  in  variable-pressure  flows. 
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